matrix_col_double.cpp